AD-A17J  814  DISTRIBUTION  ANAL VS  IS  OF  MAJOR  AND  TRACE  ELEMENTS  171 

THROUGH  SEMICONDUCTOR  L  iU)  AEROSPACE  CORP  EL  SEGUNDO 
CA  MATERIALS  SCIENCES  LAS  A  A  GALUSKA  ET  AL  15  AUG  86 
UNCLASSIFIED  TR-0886 1 6935-85 i -1  SD-TR-86-52  F7G  28718  NL 


'CROCOPy  RESOLUTION  TEST  CHART 

MTIONII  bureau  of  STANDARDS. j%3-a 


REPORT  SD-TR-86-53 


Distribution  Analysis  of  Major  and  Trace  Elements 
Through  Semiconductor  Layers  of  Changing 
Matrix  Using  SIMS 


t 

D 

> 


A.  A.  GALUSKA 
Materials  Sciences  Laboratory 
Laboratory  Operations 
The  Aerospace  Corporation 
El  Segundo,  CA  90245 

and 

G.  H.  MORRISON 
Department  of  Chemistry 
Cornell  University 
Ithaca,  NY  14853 


15  August  1986 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


Prepared  lor 

SPACE  DIVISION 
AIR  IORC  I  SYS  I  IMS  COMMAND 
I  os  Angeles  Air  force  Station 
P.O.  Box  92960.  Worldway  Postal  Ccntci 
1  os  Anucles.  CA  90009  2960 


This  report  was  submitted  by  The  Aerospace  Corporation,  El  Segundo,  CA 
90245,  under  Contract  No.  F04701-85-C-0086  with  the  Space  Division,  P. 0.  Box 
92960,  Worldway  Postal  Center,  Los  Angeles,  CA  90009-2960.  It  was  reviewed 
and  approved  for  The  Aerospace  Corporation  by  R.  W  .  Fillers,  Director, 
Materials  Sciences  Laboratory.  Lt  Ethel  Dotts,  SD/YUE,  was  the  project 
officer  for  the  Mission-Oriented  Investigation  and  Experimentation  (MOIE) 
Program. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PAS)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it 
will  be  available  to  the  general  public,  including  foreign  nationals. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the 
report's  findings  or  conclusions.  It  is  published  only  for  the  exchange  and 
stimulation  of  ideas. 


o ?l  J „/  8(. 


ETHEL  DOTTS,  Lt,  USAF 
MOIE  Project  Officer 
SD/YUE 


JOSEPH  HESS,  GM-15 

Director,  AFSTC  West  Coast  Office 

AFSTC/WCO  0L-AB 


f’t 


s 


si 


a 


$ 

'  # 
A* 


i 


UNCLASSIFIED _ 

st  cum  tv  classification  of  this  page  cpb—  Data  snt*f*d) _ 

REPORT  DOCUMENTATION  PAGE  befo^co^pletS^form 

I.  report  NUM8CR  1 2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 

SD-TR-86-53  tkitL 


|  4.  TITLE  f«iW  Sub  till*) 


TYPE  of  REPORT  *  PERIOD  COVERED 


DISTRIBUTION  ANALYSIS  OF  MAJOR  AND 
TRACE  ELEMENTS  THROUGH  SEMICONDUCTOR 
LAYERS  OF  CHANGING  MATRIX  USING  SIMS 


«.  PERFORMING  ORO.  REPORT  NUMBER 

TR-0086(6935-05)-‘ 


17.  AUTHORf  •> 


t.  CONTRACT  OR  GRANT  NUMBER^*) 


Alan  A.  Galuska  and  G.  H.  Morrison 


F04701-85-C-0086 


3.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  b  WORK  UNIT  NUMBERS 


The  Aerospace  Corporation 
El  Segundo,  CA  90245 


"Spca°^*filV2ifoFnnCE  NAM£  AND  ADOflESS 
Los  Angeles  Air  Force  Station 


IZ.  REPORT  DATE 

15  August  1986 


Los  Angeles,  CA  90009-2960 


13.  NUMBER  OF  PAGES 

46 


14.  MONITORING  AGENCY  name  «  AOORESSfi/  diltarant  /ram  Controlling  Ot Hem)  I  IS.  SECURITY  CLASS,  (at  thla  rapart) 


Unclassified 


IS  A.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


IB.  DISTRIBUTION  STATEMENT  (at  ml*  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (at  Ota  abatraet  antarad  In  Btook  30,  It  (Uttar ant  tram  Raport) 


is.  supplementary  notes 


I  is.  KEY  WOROS  (Cantlnua  on  rararaa  at da  If  naeaaaary  ltd  tdanttty  by  black  tnmibar) 


Secondary  Ion  Mass  Spectrometry 

Thin  Film  and  Interface  Characterization 


\  I  20.  ABSTRACT  (Contlnuo  on  rovoroo  oidm  it  nmemoomry  mnd  Identity  by  block  mtmbor) 


^Secondary  ion  mass  spectrometry  (SIMS)  is  often  used  to  monitor  elemental 
distributions  in  solids  and  at  solid  interfaces^  The  technique  is  highly 
sensitive  for  most  elements  and  has  good  depth  resolution.  However,  the 
complexity  of  the  sputtering  event  has  made  quantitative  analysis  difficult, 
especially  in  multimatrix  samples  where  changing  matrix  effects  are  en¬ 
countered.  Because  of  the  increasing  importance  of  raultilayer-multimatl'i*^ 


00  F0RM  U73 

(FA  CSIMI  LCl 


I  UNCLASSIFIED _ 

\  SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWiMt  Data  Sntarad) 


<9.  KEY  WORQSrCoaUiMM*) 


CONTENTS 


I.  INTRODUCTION .  7 

II.  SIMS  QUANTIFICATION .  11 

III.  MATRIX  EFFECTS . 15 

A.  Sputtering  Yield  Approach .  15 

B.  Compositional  Approach.. . 16 

C.  Oxide  Bond  Approach... . 17 

'  IV.  MATRIX  EFFECT  CALIBRATION  FOR  A1  Ga ,  As 

MATRICES .  19 

V.  QUANTITATIVE  ANALYSIS  OF  Al^Ga^As 

|  MULTILAYERS .  25 

VI.  EXTENDED  MATRIX  EFFECT  CALIBRATION .  35 

VII.  FUTURE  TRENDS .  *7 

REFERENCES .  49 


FIGURES 


1.  A  Hypothetical  Al^Ga^^As  Superlattice . . .  8 

2.  SIMS  Depth  Profile  of  a  250  keV  ^B+  Implant  into  a 

GaAs/AlQ  i2^a0  gg^s/GaAs  Sample . . .  27 

3.  SIMS  and  RBS  Depth  Profiles  of  an  AlxGaj_xAs 

Superlattice . . . . .  29 

4.  An  Uncorrected  SIMS  Depth  Profile  of  7^As+  and 

27A12+  in  an  A^Ga^^As  Superlattice .  31 

5.  A  Matrix  Corrected  SIMS  Depth  Profile  of  a  250  keV 


11B+  Implant  into  a  GsAs/AIq  3QGaQ. 7qAs  Sample 

with  a  Matrix-Dopant  Mismatch  at  the  Interface  Located 


at  a  Depth  of  Approximately  0.69  ym  [B  cone. 

(2  x  10  ®  atoms/ cm2  full  scale)  and  A1  cone. 

(lx  1022  atoms/cm'*  full  scale)].. . . .  32 

6.  SIMS  Depth  Profile  of  a  AlxGaj_xAs  Superlattice  Doped 

with  Be  and  Si . . .  34 

7.  Influence  of  Matrix  Effects  on  the  T  Value  of  2®Si .  36 

8.  SIMS  Analysis  of  llB+,24Mg+,  and  7U+ 

Implants  Through  a  GsAs/AIq^  ^Gsq  7 As  Structure .  41 

9.  SIMS  Analysis  of  «a  Be  Plateau  at  a  Graded 

GaAs/AlQ' 2GaQ ' jAs  Interface..... . . . . .  42-45 


TABLES 


1.  Practical  Ion  Yield  Calibration  for  Al^Gaj^As 

Matrices .  21 

2.  Relative  Sensitivity  Factor  Calibrations  far 

AlxGa^_xAs  Matrices . 22 

3.  Relative  Ion  Yield  Calibration  for  AlxGaj_xAs 

Matrices .  23 

4.  Relative  Sputtering  Yield  Calibrations  For 

AlxGa i_xAs  matrices .  24 

5.  Point-by-Point  Comparison  of  A1  Concentrations 

Determined  by  SIMS  and  RBS  Analyses .  30 

6.  Linearity  of  the  Sputterlig  Yield  and  Compositional 
Approaches  to  Matrix  Calibration  for  Group  IIIA-VA 

Compound  Matrices  under  O2  Bombardment . 37 

7.  Dependence  of  Rt  Versus  Bond  Energy  Line 

Slopes  on  First  Ionization  Potential .  39 


I.  INTRODUCTION 


In  Che  past,  semiconductor  technology  was  largely  limited  to  silicon 
technology.  Electronic  devices  were  rather  large  and  simple  in  construc¬ 
tion.  Current  trends  in  semiconductor  technology  involve  the  fabrication  of 
smaller  (<1.0  Pm)  and  more  complex  electronic  structures  composed  of  an 
increasing  variety  of  materials.  In  particular,  microwave  and  optoelectronic 
devices  composed  of  group  IIIA-VA  and  IIB-VIIA  compound  semiconductors  are 
rapidly  being  developed  using  growth  techniques  such  as  molecular  beam  epitaxy 
(M3E)  and  metal-organic  chemical  vapor  deposition  (MOCVD).  An  increasing 
number  of  these  compound  semiconductor  devices  are  making  use  of  superlattice 
structures,  alternating  thin  lattice-matched  semiconductor  layers  of  varying 
composition  (Fig.  1). 

There  is  a  tremendous  need  for  major  and  trace  element  quantification 
through  these  semiconductor  structures.  There  are  a  few  surface  analysis 
techniques,  such  as  Rutherford  backs catte ring  spectrometry  (RBS)  and  Auger 
electron  spectrometry  (AES),  which  have  the  necessary  depth  resolu¬ 
tion  (<  100  A)  and  sensitivity  to  examine  the  distribution  of  the  major  ele¬ 
ments.  However,  secondary  ion  mass  spectrometry  (SIMS)  is  the  only  technique 
capable  of  simultaneously  measuring  major  and  also  trace  (~  1  ppm)  elemental 
depth  distributions.  Consequently,  SIMS  has  great  potential  for  charac¬ 
terization  of  these  semiconductor  structures. 

Because  of  the  complexity  of  the  sputtering  and  ion  emission  processes, 
quantitative  SIMS  analyses  are  difficult  to  obtain.  Thus  the  full  potential 
of  SIMS  has  not  been  realized.  Every  measurement  is  influenced  by  the 
uncertainty  resulting  from  statistical  fluctuations  in  the  measured  secondary 
ion  current  and  systematic  uncertainties  related  to  experimental  conditions  as 
well  as  sample  characteristics.  For  trace  elements  in  homogeneous  matrices, 
these  sources  of  error  can  generally  be  eliminated  or  reduced  so  that 
quantitative  measurements  (~  15 Z)  can  be  made.  Conversely,  quantitative 
analysis  in  heterogeneous  matrices  has  been  severely  limited  because  of  the 
uncharacterized  variation  of  secondary  ion  yields  and  sputtering  yields  with 
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Fig  1. 


A  Hypothetical  Al^Ga^^As  Superlattice.  The  thickness  Z  of  the 
layers  can  vary  from  several  angstroms  to  several  microns  while  x 
can  be  varied  from  0  to  l.®® 
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sample  matrix  (matrix  effects).  A  depth  profile  of  a  sample  composed  of  more 
than  one  matrix  (in  particular,  layers  of  different  matrices  stacked  on  top  of 
each  other)  becomes  distorted  by  the  variation  of  secondary  ion  yields  and 
sputtering  yields  with  matrix.  Signal  and  time  are  no  longer  proportional  to 
concentration  and  depth,  respectively. 

Despite  these  difficulties,  the  proper  calibration  of  ion  yields  and 
sputtering  yields  with  matrix  can  make  quantitative  SIMS  depth  profiling 
possible  for  these  multilayer-multimatrix  samples.  In  this  report,  procedures 
for  quantifying  SIMS  depth  profiles  through  layers  of  IIIA-VA  semiconductors 
will  be  described.  The  basics  of  SIMS  quantification  and  theories  describing 
matrix  effects  will  be  reviewed,  matrix  effect  calibrations  and  depth  profile 
corrections  for  AlxGa^_xAs  multilayers  will  be  examined,  a  method  of  extending 
these  calibrations  to  the  remaining  IIIA-VA  semiconductors  will  be  presented, 
and  future  trends  in  matrix  effect  calibration  and  depth  profile  correction 
will  be  discussed. 


II.  SIMS  QUANTIFICATION 


When  a  solid  surface  is  bombarded  by  energetic  (5-20  keV)  primary  ions, 
the  elements  present  are  sputtered  primarily  as  neutrals,  although  positive 
ions,  negative  ions,  and  polyatomic  ions  are  also  formed.  Typically,  only  the 
secondary  ion  intensities  of  the  monomeric  ions,  either  M+  or  M~,  are  used  in 
SIMS  elemental  quantification.  The  relationship  between  the  number  of  analyte 
ions  detected,  1^  (counts),  and  the  isotopically  corrected  analyte 
concentration,  Cv  (atoras/cm^),  in  the  sample  surface  can  be  defined  as 

ra ' 'm  cm  sm  j  \>  T  S'1  (1) 

where  is  the  useful  ion  yield,  SM  is  the  sputtering  yield,  J  is  the 
current  density  (ions/cm^-sec) ,  and  A  is  the  analysis  area  (cm^).  T  is  the 
analysis  time  (sec),  and  N  is  the  atomic  density  (atoms/cm^)  of  the  sample 
matrix.  The  useful  ion  yield  is  defined  as  the  ratio  of  the  number  of  analyte 
ions,  M*,  detected  to  the  number  of  analyte  atoms  sputtered  from  the 
specimen.  Useful  ion  yields  are  also  defined  as  the  product  of  the 
probability  of  M*  formation  (8*)  and  the  transmission  efficiency  (n ) .  The 
sputtering  yield  is  defined  as  the  average  number  of  atoms  sputtered  by  each 
i  -dent  primary  ion.  It  is  assumed  that  the  elemental  components  of  the 
sputtering  yield  are  proportional  to  the  bulk  elemental  concentrations. 

The  difficulty  in  solving  Eq.  (1)  for  analyte  concentration  is  that 
useful  ion  yields  and  sputtering  yields  are  influenced  by  a  variety  of 
factors,  and  can  vary  drastically  from  sample  to  sample.  Useful  ion  yields 
can  vary  by  orders  of  magnitude  for  different  elements  in  the  same  matrix  and 
for  the  same  element  in  different  matrices. In  addition,  useful  positive 
ion  yields  are  increased  by  the  presence  of  electronegative  elements  such  as 
oxygen,'’  and  useful  negative  ion  yields  are  increased  by  the  presence  of 
electropositive  elements  such  as  cesium.^  Sputtering  yields  are  also 
different  for  different  elements  in  the  same  or  different  sample  matrices. 
Sputtering  yields  are  influenced  by  the  following:  primary  ion  mass,  energy 


and  density,  mass  of  the  target  atoms,  sample  crystallinity  and  orientation, 
surface-binding  energy,  and  the  angle  of  incidence  between  the  primary  ion 
beam  and  the  sample  surface.  Consequently,  useful  ion  yields  and 
sputtering  yields  must  be  determined  for  every  element  and  sample  matrix.  In 
addition,  because  of  the  fluctuation  of  experimental  conditions  with  time, 
these  values  must  be  determined  each  time  quantification  is  performed  in  the 
same  sample  matrix. 

For  pure  elemental  matrices,  sputtering  yields  can  be  determined  experi¬ 
mentally  from  the  sputtering  rate  z  (cm/sec).  After  sputtering,  the  depth  of 
the  sputtered  crater  is  measured  using  a  profilometer  or  an  optical  inter¬ 
ferometer.  If  this  depth  is  divided  by  the  sputtering  time  interval,  the 
average  sputtering  rate  can  be  obtained.  If  the  approximation  is  made  that 
the  crater  walls  behave  as  step  functions,  the  sputtering  yield  may  be 
determined  from 


S  -  z  A  N/i  (2) 

P 

where  A  is  the  crater  area  (cm^),  and  ip  is  the  primary  beam  current 
( ions/sec) . 


In  matrices  of  more  than  one  element,  differences  in  elemental  sputtering 
yields  (preferential  sputtering)  result  in  composition  changes  in  the  near 
surface  region.*®'**  Elements  which  sputter  slowly  are  enriched  at  the 
surface  while  elements  which  sputter  rapidly  are  depleted  from  the  surface. 


These  compositional  changes  usually  occur  over  a  short  depth  increment 
(approximately  twice  the  range  of  the  primary  ion)  after  which  a  steady  state 
condition  is  reached.  At  this  point,  the  elemental  components  of  the 
sputtering  yield  are  proportional  to  the  bulk  elemental  concentrations.  Once 
this  steady  state  condition  has  been  reached,  Eq.  (2)  can  be  incorporated  into 
Eq.  (1)  to  obtain  a  more  practical  form  of  the  ion  intensity  equation. 
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Equations  1  and  3  can  now  be  solved  for  the  analyte  concentration  pro¬ 
viding  t  can  be  accurately  determined.  Useful  ion  yields  can  be  approximated 
M 

using  semitheoretical  techniques 9 12-16  ^ut  they  are  most  accurately  determined 
using  homogeneous  bulk  standards  or  ion  implanted  standards^-^  fabricated 
from  the  same  matrix  as  the  unknown  sample. 

Although  practical  ion  yields  can  be  determined  from  standards,  they  can 
change  quite  dramatically  between  analyses.  Secondary  ion  collection  and 
instrumental  transmission  are  influenced  by  operator  adjustment  and  inde¬ 
terminate  instrumental  fluctuations  that  occur  during  routine  operation.  To 
reduce  this  effect,  it  is  common  practice  to  use  relative  sensitivity  factors 
(RSF),  as  given  by  Eq.  (4),  in  which  the  t  of  an  analyte  is  ratioed  to  that  of 
a  reference  element  of  known  concentration  in  the  sample  matrix. 


The  normalization  procedure  is  designed  to  minimize  the  influence  of  instru¬ 
mental  variations. 


III.  MATRIX  EFFECTS 


The  fabrication  of  external  and  internal  standards  by  ion  implantation 
has  provided  an  accurate  (15%  relative  standard  deviation)  means  of 
quantifying  trace  element  distributions  in  homogeneous  matrices.  However, 
because  of  the  uncharacterized  variation  of  secondary  ion  yields  and 
sputtering  yields  with  matrix  composition  (matrix  effects),  the  quantitative 
SIMS  analysis  of  heterogeneous  matrices  remains  a  problem.  To  overcome  this 
difficulty,  the  mechanisms  underlying  matrix  effects  must  be  more  clearly 
understood.  Several  factors  are  known  to  influence  matrix  effects.  These 
include  the  bulk  matrix  composition  and  the  surface  concentration  of  reactive 
species  (electronegative  elements  such  as  oxygen  or  electropositive  elements 
such  as  cesium),  which  can  either  be  adsorbed  from  the  sample  chamber  during 
analysis  or  implanted  as  a  primary  ion.  Several  explanations  for  matrix 
effects,  which  depend  to  varying  extents  on  the  bulk  composition  and  the 
reactive  species  concentrations  have  been  proposed.  The  sputtering  yield, 
compositional,  and  oxide  bond  approaches  to  matrix  effect  calibration  will  be 
presented  and  briefly  discussed. 

A.  SPUTTERING  YIELD  APPROACH 

The  sputtering  yield  approach  to  matrix  effect  calibration  is  based  on 
the  fact  that  positive  and  negative  ion  yields  are  greatly  enhanced  by  the 
incorporation  of  oxygen  or  cesium,  respectively,  in  the  sputtered  surface. 
Deline  et  al.  J  state  that  for  reactive  beam  sputtering,  "ion  yield  matrix 
effects  arise  solely  from  variations  in  the  substrate  sputtering  yields." 

They  justify  this  statement  by  asserting  that  implanted  oxygen  and  cesium  are 
responsible  for  the  observed  matrix  effects,  and  that  the  equilibrium  concen¬ 
tration,  [R] ,  of  these  reactive  species  is  determined  by  the  sputtering  yield, 
S. 

The  relationship  between  [ R ]  and  S  is  derived  from  equations  describing 

24  25 

ion  implantation  during  sputtering:  * 


(6) 


C  -  [R]/( [R]  +  N)  -  a/(a  +  S) 
s 

where  a  is  Che  acconmodation  coefficient  (which  equals  1  —  B,  where  B  is  Che 
backscattering  coefficient);  N  is  Che  atomic  density  of  the  sample  surface; 
and  Cs  is  the  atomic  fraction  of  implanted  projectiles  in  the  surface.  In 
their  declaration.  Deline  et  al.  assume  that  a  -  1,  that  N  and  S  are  constant, 
and  that  there  is  no  preferential  sputtering.  In  their  experiments,  they  make 
the  further  assumption  that  when  primary  current  density  is  constant,  the 
sputtering  rate  is  approximately  proportional  to  the  sputtering  yield 
regardless  of  N. 

The  result  of  these  approximations  is  the  relationship 


t  -  [R]«(l/S)y  -  (1/S)y  _ 


const 


where  y  is  a  constant  determined  experimentally  for  each  analyte.  Plots  of 
log  (t)  versus  log  (1/z)  obtained  by  Deline  et  al. 21-23  and  others26  for 
various  elements  have  been  roughly  linear.  However,  the  sputtering  yield 
approach  has  been  severely  criticized2^-30  because  of  the  many  unjustified 
assumptions  that  are  made.  Various  attempts  to  experimentally  determine  the 
relationship  between  [R]  and  S  using  Auger  spectroscopy  have  both  supported 
and  contradicted  the  proposed  relationship. It  i3  unclear  when  the 
sputtering  yield  approach  applies  or  does  not  apply. 

B.  COMPOSITIONAL  APPROACH 

The  compositional  approach  is  based  on  experiments,  performed  on  selected 

ternary  metal  alloys  and  silicates,  which  indicate  a  linear  relationship 

31—33 

between  practical  ion  yields  and  matrix  composition.  Using  oxygen 

flooding  in  conjunction  with  a  low  current  density  Ar+  primary  beam,  the 
sample  surfaces  were  saturated  with  oxygen.  Thus,  the  effects  of  reactive 
species  enhancement  were  standardized.  The  remaining  matrix-induced 
variations  of  practical  ton  yields  were  then  explained  using  the  expression 
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(8) 


M 


v n  p  c 
i-1  M,i  i 


where  is  the  practical  ion  yield  of  element  M;  ^  is  a  parameter  repre¬ 
senting  the  influence  of  element  i  on  the  ion  yield  of  element  M;  and  is 
the  atomic  fraction  of  element  i. 

For  quantitative  analysis  based  on  this  approach,  relative  sensitivity 
factors  were  obtained  using  Eq.  (8)  to  described  the  practical  ion  yields  of 
both  the  analyte  and  the  reference.  Relative  sensitivity  factors  obtained  by 
such  a  calculation  were  shown  to  agree  with  the  experimentally  determined 
values  to  an  accuracy  of  1  to  72  RSD. ^  Despite  the  success  of  this  method  in 
the  cases  studied,  it  has  not  been  successfully  applied  to  other  types  of 
samples.  Although  the  reason  for  the  limited  use  of  the  method  has  not  been 
stated  in  the  literature,  it  is  probably  related  to  the  difficulty  of 
saturating  nonmetallic  surfaces  with  oxygen. 

C.  OXIDE  BOND  APPROACH 

The  oxide  bond  explanation  of  matrix  effects  is  based  on  the  work  of  Yu 
and  Reuter. They  performed  SIMS  analyses  on  various  binary  metal  alloys 
under  both  0*  bombardment  and  Ar+  bombardment  with  oxygen  flooding.  After  the 
SIMS  analyses,  the  alloys  were  then  analyzed  with  X-ray  photoelectron  spec¬ 
troscopy  (XPS)  to  determine  the  types  of  metal  oxides  formed  and  the  extent  of 
oxidation. 


The  XPS  data  from  the  pure  elemental  metals  revealed  that  oxidation  was 
more  extensive  for  those  metals  that  form  metal  oxides  with  large  negative 
free  energies  of  formation  than  those  metals  that  form  oxides  with  less 
favorable  free  energies  of  formation.  The  addition  of  the  second  alloy  com¬ 
ponent  resulted  in  significant  changes  in  oxidation,  ionization  probability, 
and  ion  energy  distribution  for  both  components.  The  observed  trends  can  be 
summarized  by  two  general  rules.  First,  for  an  alloy  A-B,  where  A  forms  a 
stronger  oxide  bond  than  B,  the  presence  of  A  enhances  the  oxidation  and  also 
the  ion  yield  of  B,  whereas  the  presence  of  B  suppresses  the  oxidation  and 
also  the  ion  yield  of  A.  Second,  the  presence  of  A  sharpens  the  energy 
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distribution  of  B  ,  whereas  the  presence  of  B  broadens  the  energy  distribution 
of  A+.  Both  the  trends  in  ion  yield  and  ion  energy  distribution  were  attri¬ 
buted  to  the  relative  extent  of  surface  oxidation.  The  influence  of  surface 
oxidation  has  been  attributed  both  to  changes  in  the  work  function  of  the 
surface  and  also  to  changes  in  the  bond  breaking  mechanism  of  the  surface. ^ 
While  these  rules  were  firmly  supported  by  data,  they  were  never  applied  to 
any  type  of  matrix  effect  calibration  for  quantitative  SIMS  analysis. 
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IV.  MATRIX  EFFECT  CALIBRATION  FOR  AlxGa1_xAs  MATRICES 

There  have  been  only  a  few  investigations  of  the  influence  of  matrix 
effects  on  the  ion  yields  and  sputtering  yields  of  IIIA-VA  semiconductors. 

For  homologous  series  of  compounds,  such  as  Aj^Ga^^As  and  GaAsxPj_x  matrices, 
under  both  0^  and  Cs  bombardment,  KatzJ  observed  that  sputtering  yields  were 
linearly  related  to  the  changing  matrix  composition.  Similarly,  for  AlxGaj_xAs 
matrices  under  0*  bombardment,  others30,38  have  observed  that  the  useful  ion 
yields  of  major  and  trace  elements  were  linearly  related  to  the  matrix 
composition  defined  by  x.  Although  these  investigations  indicated  that  linear 
relationships  existed,  precise  and  reproducible  matrix  effect  calibration 
lines  could  not  be  obtained  because  of  the  indeterminate  variation  of  instru¬ 
mental  parameters  that  alter  the  observed  ion  yields  and  sputtering  yields 
between  analyses.  A  method  of  normalizing  these  instrumental  parameters  was 
required. 

The  significance  of  instrumental  variations  and  procedures  for  nor- 
malizing  their  influence  were  examined  by  Galuska  and  Morrison.  Using 
an  0*  primary  ion  beam,  they  performed  a  series  of  analyses  on  groups  of  ion 
implanted  AlxGaj  As  matrices.  To  ensure  nearly  identical  analysis  conditions 
during  each  analysis  period,  groups  of  samples,  including  in  each  case  GaAs , 
were  inserted  simultaneously.  These  samples  were  sequentially  analyzed 
without  manipulating  any  instrumental  parameters.  While  minimizing  instru¬ 
mental  variations  in  this  way,  the  useful  ion  yields  (t)  of  the  trace  and 
major  elements  were  obtained.  These  useful  ion  yields  were  then  ratioed  to 
produce  relative  sensitivity  factors  (RSFs)  and  relative  ion  yields  (Rts). 

Rt  -  t  /t  (9) 

x  0 

where  t  is  the  useful  ion  yield  of  the  analyte  in  a  particular  Al^Gaj^As 

matrix,  and  r  is  the  useful  ion  yield  of  the  analyte  in  GaAs  when  both 
o 

measurements  are  performed  under  nearly  identical  conditions. 


The  linearity,  precision,  and  reproducibility  of  matrix  effect  calibra¬ 
tion  lines  formed  using  t's,  RSFs  and  Rts  were  then  compared.  The  results  for 
Be,  Si,  P,  B,  and  As  are  presented  in  Tables  1,  2,  and  3.  The  calibration 
lines  obtained  using  relative  ion  yields  were  the  most  linear,  precise,  and 
reproducible.  As  expected,  the  calibration  lines  obtained  using  useful  ion 
yields  were  linear  but  exhibited  poor  precision  and  reproducibility.  Despite 
good  reproducibility  (~  24.4%)  for  a  given  data  point,  the  relative 
sensitivity  factors  resulted  in  calibration  lines  with  poor  linearity, 
precision,  and  reproducibility.  The  superiority  of  Rts  over  RSFs  for  matrix 
effect  calibration  is  readily  justified.  An  Rr  is  directly  proportional  to 
the  ion  yield  of  the  analyte ,  whereas  an  RSF  is  related  to  the  ion  yields  of 
both  the  analyte  and  the  reference.  Because  the  influence  of  instrumental 
variations  on  the  ion  yields  of  the  analyte  and  reference  can  vary 
independently,  RSFs  are  altered  by  instrumental  changes  to  a  greater  extent 
than  Rrs.  Moreover,  Rts  are  designed  to  normalize  matrix  effects  to  a 
standard  matrix,  whereas  RSFs  are  not  so  designed. 

In  addition  to  t's,  RSFs,  and  Rrs,  absolute  sputtering  yields^  and 
relative  sputtering  yields  (RSs)  were  also  determined  as  a  function  of  matrix 
composition.  As  shown  in  Table  4,  the  calibration  lines  obtained  using 
relative  sputtering  yields  (where  RS  -  SX/S(J)  were  quite  precise  and  repro¬ 
ducible,  unlike  the  absolute  sputtering  yields.  The  use  of  Rrs  and  RSs  has 
permitted  the  influence  of  matrix  changes  on  ion  yields  and  sputtering  yields 
to  be  precisely  calibrated  (~  10%  RSD)  for  the  first  time. 
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Table  1.  Practical  Ion  Yield  Calibrations  for  AlxGa^_xAs  Matrices38 


Analyte 

Week  No. 

X  Values 

Intercept 
( x  10‘7) 

Slope 
(x  10“5) 

Linear 

Correlation 

(r2) 

RSD 

Slope 

(%) 

9Be 

1 

0,0.12,0.26,0.37 

1.23 

41.5 

0.927 

32.8 

2 

o 

o 

• 

►— * 

u> 

o 

• 

ro 

h-* 

o 

• 

UJ 

116 

67.3 

0.915 

68.6 

3 

0,0.13,0.18,0.21 

61.4 

32.3 

0.995 

9.5 

28Si 

1 

36.4 

22.2 

0.989 

18.5 

2 

25.3 

15.2 

0.987 

17.0 

3 

38.9 

17.7 

0.984 

10.8 

31p 

16 

0,0.18,0.26,0.31 

0.822 

0.291 

0.936 

18.6 

32 

0,0.13,0.18,0.37 

1.17 

0.258 

1.000 

2.0 

UB 

16 

25.8 

7.05 

0.734 

53.6 

32 

28.2 

17.0 

0.994 

8.7 

75As 

1 

1.78 

0.130 

0.966 

14.6 

2 

0.932 

0.0891 

0.915 

21.0 

3 

0.940 

0.108 

0.984 

11.0 

32 

0.918 

0.0532 

0.932 

11.2 

Av  RSD  per 

Point  ■  18. 

.42 

9Be 

total 

31.4 

44.7 

0.868 

12.3 

28Si 

total 

11.0 

19.0 

0.994 

5.8 

31P 

total 

1.00 

1.00 

0.278 

0.928 

12.4 

llB 

total 

27.0 

15.3 

0.825 

18.8 

75As 

total 

0.642 

0.111 

0.722 

19.2 

Av  RSD  per  Point  ■  60. 42 


Table  2.  Relative  Sensitivity  Factor  Calibrations  for 
Al^Gai^As  Matrices™ 


Analyte 

Week  No. 

X  Values 

Intercept 
(x  10-7) 

Slope 
(x  10"5) 

Linear 

Correlation 

(r2) 

RSD 

Slope 

(Z) 

9Be/75As 

1 

0,0.12,0.26,0.37 

47.0 

890 

0.731 

33.6 

2 

0,0.13,0.21,0.37 

84.8 

573 

0.894 

23.0 

3 

0,0.13,0.18,0.21 

66.5 

950 

0.995 

4.7 

2hi/7^ 

1 

39.8 

271 

0.708 

38.2 

2 

39.4 

382 

0.903 

24.2 

3 

36.0 

711 

0.936 

37.8 

31P/75As 

16 

0,0.18,0.26,0.31 

1.58 

9.67 

0.715 

40.3 

32 

0,0.13,0.18,0.37 

1.88 

4.91 

0.827 

58.9 

HB/75As 

16 

30.7 

454 

0.945 

26.1 

32 

60.9 

509 

0.940 

17.1 

Av  RSD  per 

Point  “  10. 

4Z 

9Be/7!)As 

total 

88.7 

740 

0.694 

16.8 

28Si/75As 

total 

42.7 

357 

0.804 

17.7 

31P/75As 

total 

2.41 

4.45 

0.464 

79.2 

UB/75As 

total 

0 

34.0 

481 

0.883 

19.2 

Av  RSD  per 

Point  •  24. 

4Z 
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Table  3.  Relative  Ion  Yield  Calibrations  for  AlxGa^_xAs  Matrices88 


Analyte 

Week  No. 

Intercept 

X  Values  (x  1CT7) 

Slope 
(x  10~5) 

Linear 

Correlation 

(r2) 

RSD 

Slope 

(7.) 

9Be 

1 

0,0.12,0.26,0.37  1.0 

52.9 

0.999 

6.4 

2 

0,0.13,0.21,0.37  1.0 

42.6 

0.996 

4.3 

3 

0,0.13,0.18,0.21  1.0 

41.7 

1.000 

14.2 

28Si 

1 

1.0 

44.5 

1.000 

9.8 

2 

1.0 

40.1 

0.975 

5.7 

3 

1.0 

45.9 

1.000 

9.6 

3 1  p 

16 

0,0.18,0.26,0.31  1.0 

21.2 

0.995 

9.6 

32 

0,0.13,0.18,0.37  1.0 

23.0 

1.000 

2.5 

“b 

16 

1.0 

46.7 

0.990 

10.7 

32 

1.0 

39.5 

1.000 

8.6 

78As 

1 

1.0 

6.00 

0.933 

15.7 

2 

1.0 

6.54 

0.991 

2.8 

3 

1.0 

5.76 

0.980 

12.4 

32 

1.0 

5.32 

.  0.925 

11.0 

Av  RSD  per  Point  ■  7. 

42 

9Be 

total 

1.0 

46.3 

0.981 

13.7 

28S1 

total 

1.0 

44.7 

0.986 

5.0 

31p 

total 

1.0 

22.1 

0.942 

4.6 

total 

1.0 

46.6 

0.961 

7.0 

75As 

total 

1.0 

5.76 

0.917 

7.0 

Av  RSD  per  Point  ■  9. 

27. 

V.  QUANTITATIVE  ANALYSIS  OF  AlxGa1_xAs  MULTILAYERS 


Through  the  use  of  Rt  and  RS  calibration  lines,  the  quantitative  SIMS 
analysis  of  AlxGa^_xAs  multilayers  (superlattices)  has  become  possible.  The 
only  requirement  is  that  a  GaAs  standard  be  analyzed  under  the  same  conditions 
as  the  sample  of  interest.  The  concentration  of  analyte  at  each  point  of  a 
depth  profile  Cp  (atora/cm^)  can  be  determined  using  Eq.  (10). 


C  -  I  /RT  t  z  A  (10) 

p  p  X  o  p  o 


where  RTx  is  the  relative  ion  yield  determined  from  a  calibration  line  for  the 
appropriate  value  of  x,  Tq  is  the  useful  ion  yield  of  the  analyte  in  the 
standard  matrix  (GaAs),  Ip  is  the  signal,  and  zp  is  the  depth  increment 
associated  with  each  point  of  the  depth  profile.  Similarly,  the  erosion  rate 
at  each  point  of  a  depth  profile  zp  can  be  determined  using  Eq.  (11). 


z  -  RS 
p  x 


o 

z 

o 


N  /N 

o  X 


(11) 


where  RSX  is  the  relative  sputtering  yield  determined  from  a  calibration  line 

for  the  appropriate  value  of  x;  z  is  the  sputtering  rate  of  the  standard 

o 

matrix  (GaAs);  and  NQ  and  are  the  atomic  densities  of  the  standard  matrix 
and  the  sample,  respectively. 

The  application  of  Eqs.  (10)  and  (11)  requires  a  knowledge  of  the  matrix 
composition  (x)  at  each  data  point.  This  information  is  not  readily  available 
from  an  uncalibrated  SIMS  depth  profile.  However,  a  program,  superlattice  and 
interface  calibration  (SLIC),^  has  been  developed  that  determines  x  at  each 
data  point  using  the  ^As+  calibration  line,  the  RS  calibration  line,  and  the 
fact  that  the  As  concentration  is  constant  in  A^Ga^^As  matrices.  If  the 
value  of  Rt^  in  Eq.  (10)  is  expressed  in  terms  of  the  equation  for  the  ^As+ 
calibration  line  (Rt.  ■  x  M.  +  1 ,  where  M.  D  is  the  slope  of  the 
line),  the  equation  can  be  manipulated  to  the  following  form: 
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X 


(12) 


{I.  /CA  *  A  *t  *z  )  -  l}/M.  _ 

As  ,p  As  ,p  o  As  ,o  p  As,Rr 


Except  for  zp,  the  values  for  all  the  variables  on  the  right  side  of  Eq.  (12) 

can  be  readily  determined.  Because  z_  is  directly  related  to  z  (z_  »  z  • 

P  P  P  P 

Tp,  where  Tp  is  the  number  of  seconds  per  point),  it  can  be  expressed  in  terms 
of  Eq.  (11).  If  RSX  is  replaced  with  the  corresponding  equation  for  the  RS 
calibration  line,  Eq.  (13)  is  obtained 


<mrs  • 5  ♦  *>  •  K  ■  • 

CMrf  •  *  *  D  •  s0  •  T 


where  is  the  slope  of  the  RS  calibration  line,  and  M^o  is  the  slope  of  the 

corresponding  relative  erosion  rate  calibration  line  (M_o  *  MR<>  •  Nx/Nq). 

Once  again  the  values  of  zp  and  x  are  the  only  unknown  quantities. 

SLIC  initially  use3  the  assumption  that  x  »  0.  The  corresponding  value 
of  zp  is  then  determined  from  Eq.  (13)  and  employed  in  Eq.  (12)  to  obtain  a 
better  approximation  of  x.  This  process  is  reiterated  until  the  values  of  x 
and  Zp  converge.  In  this  manner,  the  matrix  composition  and  depth  at  each 
point  of  a  depth  profile  are  determined  through  matrix  gradients,  interfaces, 
and  plateaus.  The  dopant  profiles  are  then  corrected  for  matrix  changes  at 
each  data  point  using  Eq.  (10)  and  the  appropriate  dopant  calibration  lines. 

A  250  keV  implant  into  a  layered  sample  (GaAs/Al^  ^ 2^.88^ 
provides  a  good  example  for  the  correction  of  matrix  effects.  In  Fig.  2a,  the 
uncorrected  SIMS  depth  profile  of  ^B+  and  ^As+  in  this  sample  is  pre¬ 
sented.  The  first  interface  occurs  at  33  time  units,  whereas  the  second 
interface  has  not  been  reached  in  this  profile.  Because  of  matrix  effects, 
the  vertical  and  horizontal  scales  are  no  longer  linearly  related  to 
concentration  and  depth,  respectively.  The  peak  in  the  ^B+  profile  at  34 
time  units  is  a  good  example  of  a  distortion  introduced  by  the  changing  matrix 
effects  at  the  interface.  In  Fig.  2b,  the  horizontal  scale  has  been 
transformed  into  depth,  and  the  concentration  of  Al  at  each  point  has  been 
calculated  from  the  ^As+  profile.  The  ^B+  profile  was  then  quantified.  An 
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A1  concentration  of  2.4  x  10^  atom/cm^  ±  6.02  was  calculated  for  the  AlxGa^_ 
xAs  layer,  which  agrees  quite  well  with  the  value  expected  from  MBE  growth 

parameters.  The  correction  procedure  has  also  removed  the  distorted  interface 
region  of  the  11B+  profile,  transforming  the  profile  into  the  slightly 
distorted  Gaussian  shape,  which  was  expected. 

The  most  crucial  factor  influencing  the  quality  of  the  matrix  corrections 
performed  by  SLIC  is  the  accuracy  with  which  the  matrix  structure  can  be 
determined.  This  accuracy  has  been  checked  with  RBS  analysis,  which  is  accu¬ 
rate  to  about  102  RSD  and  has  a  detection  limit  of  12  for  AlxGa^_xAs 
matrices.  In  Fig.  3,  a  complex  A^Ga^^As  sample  was  analyzed  by  SIMS  and 
RBS,  respectively.  The  A1  concentration  at  specific  regions  of  the  super¬ 
lattice  was  determined  using  SLIC  and  standard  RBS  techniques.^®  As  shown  in 
Table  5,  the  values  determined  by  the  two  techniques  agree  quite  well.  In 
fact,  the  two  sets  of  data  are  not  statistically  different  at  a  952  confidence 
level. 

In  addition  to  matrix  composition,  the  point-by-point  correlation  between 
matrix  structure  and  dopant  distribution  is  critical.  Small  differences 
between  the  actual  and  the  measured  matrix  structure  can  significantly  influ¬ 
ence  the  corrections  performed  on  dopant  distributions.  As  apparent  in 
Fig.  4,  the  ^As+  signal  tracks  the  A1  distribution  quite  well.  This  corre¬ 
lation  is  used  by  SLIC  to  precisely  determine  matrix  structures.  However, 
errors  can  still  result  from  the  correlation  between  matrix  structure  and 
dopant  distribution.  For  example,  there  is  a  dead  time  between  the  measure¬ 
ment  of  the  matrix  signal  and  the  dopant  signal  for  each  point  of  a  depth 
profile.  When  this  dead  time  is  large  compared  to  the  abruptness  of  the 
matrix  changes,  the  matrix  and  dopant  signals  will  be  obtained  from  different 
matrix  regions.  Consequently,  the  dopant  signal  will  be  calibrated  for  the 
wrong  matrix  composition.  Such  an  error  Is  shown  in  Fig.  5.  Fortunately, 
this  type  of  error  can  usually  be  avoided  by  minimizing  the  dead  time  and 
reducing  the  sputtering  rate. 

The  full  utility  of  SLIC  can  be  appreciated  when  very  complex  samples  are 
analyzed.  An  uncorrected  SIMS  depth  profile  of  a  complex  AlxGa^_xAs  super- 
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TIME  (arbitrary  units) 

(b) 

ENERGY  (MeV) 


Table  5.  Point-by-Point  Comparison  of  A1  Concentrations 
Determined  by  SIMS  and  RBS  Analyses^’ 


Point 

A1  Concentration  (X 

RBS 

10^*  atom/cm^) 

SIMS 

Deviation 

1 

12.0 

13.0 

-1.0 

2 

4.2 

3.9 

+0.3 

3 

5.1 

4.8 

+0.3 

4 

5.1 

4.7 

+0.4 

5 

5.2 

4.7 

+0.5 

6 

11.0 

11.0 

0.0 

7 

1.0 

1.3 

-0.3 

8 

4.9 

5.7 

-0.8 

AVG  Deviation  ■  -0.075 


~0  50  100 

TIME  (arbitrary  units) 

An  Uncorrected  SIMS  Depth  Profile  of  ^As+  (~  — )  and  ^A1 
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an  Al  Ga1_xAs  Superlattice 
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lattice  grown  by  MBE  with  Be  and  Si  dopants  is  presented  in  Fig.  6a.  In  Fig. 
6b ,  SLIC  has  been  used  to  quantify  the  elemental  distributions  in  this 
sample.  In  the  uncorrected  profile,  both  the  Be  and  Si  distributions  follow 
the  ^As+  signal  because  of  the  changing  matrix  effects.  Upon  calibration, 
both  the  Be  and  Si  distributions  have  changed  substantially.  As  expected  from 
the  MBE  growth  conditions,  the  Be  concentration  generally  increases  as  the  Al 
concentration  decreases  and  vice  versa.  In  addition,  excluding  the  surface 
buildup,  the  Si  distribution  has  generally  leveled  out  at  4  x  10  atom/cm  . 
Without  SLIC,  this  type  of  analysis  could  not  have  been  made.  At  the  present, 
SLIC  represents  the  only  successful  attempt  at  using  SIMS  to  quantify  ele¬ 
mental  distributions  through  multilayers  of  changing  matrix. 
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VI.  EXTENDED  MATRIX  EFFECT  CALIBRATION 


Although  the  initial  attempts  at  matrix  effect  calibration  and  correction 
have  been  successful  for  AlxGa^_xAs  matrices,  the  general  application  of  such 
calibration  procedures  to  a  greater  variety  of  matrices  requires  a  better 
understanding  of  the  fundamental  processes  involved.  During  their  attempts  to 
extend  the  AlxGaj_xAs  matrix  effect  calibration  lines  to  other  IIIA-VA  semi¬ 
conductors  (GaQ^jIn^ 53AS ,  GaSb,  InSb,  InP,  Gap,  and  AlxGaj_xAs  matrices  were 
analyzed),  Galuska  and  Morrison^ *  examined  the  applicability  of  the  sputtering 
yield,  compositional,  and  oxide  bond  explanations  of  matrix  effects.  As  in 
previous  analyses,  an  0+  primary  ion  beam  was  used.  To  reduce  experimental 
errors,  practical  ion  yields  and  sputtering  yields  were  replaced  with  Rr  and 
RS,  respectively. 

The  sputtering  yield  hypothesis  was  evaluated  directly  using  Eq.  (7). 

The  compositional  approach  was  examined  in  conjunction  with  the  oxide  bond 
approach.  The  parameter  PM^  from  Eq.  (8)  was  treated  as  a  measure  of  the 
affinity  of  element  i  for  oxygen.  The  values  of  Pjj  ^  were  approximated  from 
the  average  bond  energies  of  the  diatomic  metal  oxides. 

According  to  the  compositional  theory,  plots  of  Rr  versus  the  average 
matrix-oxygen  bond  energy  should  yield  straight  lines.  Similarly,  according 
to  the  sputtering  yield  hypothesis,  plots  of  log  (Rx)  versus  log  (1/RS)  should 
also  yield  straight  lines.  The  experimental  data  for  ^Si+,  plotted  using  the 
two  different  approaches,  are  {Presented  in  Figs.  7a  and  7b.  As  apparent  in 
these  figures,  the  compositional  theory  satisfies  the  linearity  criteria  much 
better  than  does  the  sputtering  yield  approach.  This  fact  is  emphasized  in 
Table  6  where  the  linear  correlations  and  standard  deviations  of  the  lines  for 
several  different  elements  are  compared.  In  each  case,  excellent  linearity 
and  precision  were  observed  when  using  the  compositional  approach,  and  poor 
linearity  and  precision  were  observed  when  using  the  sputtering  yield 
approach.  This  linearity  persisted  whether  the  analyses  were  performed 
under  0*  bombardment  or  Ar+  bombardment  with  oxygen  flooding.  These 
experiments  demonstrate  that  the  affinity  of  a  matrix  for  oxygen  seems  to  be 


LOG  (1 /relative  sputtering  yield) 

(b) 


MATRIX  BOND  ENERGY  TO  OXYGEN  (kcal/mole) 

Fig.  7.-  Influence  of  Matrix  Effects  on  the  t  Value  of  28S1.  (a)  Sputtering 

yield  approach;  (b)  compositional  approach  (Al„Ga ,  As-O,  InP-  a 
InSb-  o  .  GaP-  0  ).41  X  i“X 


Che  critical  facCor  (at  least  for  IIIA-VA  semiconductors  under  0*  bombardment) 
determining  the  variation  of  ion  yields  with  matrix.  This  oxygen  affinity  is 
determined  by  Che  matrix  composition.  Although  the  sputtering  yield  does 

Influence  the  amount  of  oxygen  available  for  bonding,  it  does  not  determine 
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the  final  oxygen  content  of  the  surface. 

To  determine  a  general  mechanism  for  matrix  effects,  it  is  Important  to 
understand  how  matrix  effects  influence  the  t's  of  one  particular  element 
versus  another.  One  can  gain  an  insight  into  this  phenomenon  by  comparing  the 
slopes  of  the  calibration  lines  to  the  first  ionization  potentials  of  the 
respective  analytes.  Such  a  comparison  is  presented  in  Table  7.  For  those 
elements  that  fall  in  the  same  column  of  the  periodic  table  [(B,  Al,  Ga) ,  (P, 
As,  Sb) ,  and  (Be  and  Mg)],  a  general  relationship  is  observed.  Elements  with 
larger  first  ionization  potentials  yield  calibration  lines  of  steeper  slope 
chan  those  with  smaller  first  ionization  potentials.  This  indicates  that 
elements  that  have  larger  ionization  potentials  (small  ionization  probability) 
are  Influenced  by  matrix  effects  (the  amount  of  bound  oxygen)  to  a  greater 
extent  chan  those  that  have  smaller  ionization  potentials  (large  ionization 
probability).  As  one  might  expect,  this  relationship  is  very  similar  to  that 
commonly  observed  for  elemental  t's  in  a  single  matrix  enhanced  by  oxygen 
bombardment  or  flooding:  the  t's  of  elements  with  small  ionization 
probabilities  are  enhanced  to  a  greater  extent  by  the  presence  of  oxygen  than 
those  of  elements  with  large  ionization  probabilities. 

The  excellent  linearity  of  the  Rr  versus  matrix-oxygen  bond  energy  cali¬ 
bration  lines  and  the  relationship  of  the  lines  to  elemental  ionization 
potentials  can  be  used  to  improve  the  quality  of  both  qualitative  and 
quantitative  SIMS  analysis.  Using  the  relationships  expressed  in  these  matrix 
calibration  lines,  one  can  anticipate  when  matrix  effects  will  be  a  problem, 
and  how  they  may  distort  depth  profiles  of  layered  multimatrix  samples.  For 
example,  a  largely  distorted  Gaussian  depth  distribution  would  be  expected  for 
the  SIMS  analysis  of  an  ^B+  implant  through  a  GaAs  layer  into  an  Al.3QGa.7gAs 
layer  because  **B+  emission  is  very  sensitive  to  the  presence  of  oxygen,  i.e., 
matrix  effects,  and  Al.3gGa.7gAs  has  a  much  greater  oxygen  affinity  than 
GaAs.  Alternatively,  a  smaller  distortion  would  be  expected  for  a  Mg 
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implant  and  an  even  smaller  distortion  for  a  'Li+  implant  into  such  a 
structure.  The  sensitivity  of  the  t's  of  these  elements  to  matrix  effects 
decreases  going  from  * *B+  to  24Mg+  to  7Li+.  The  SIMS  analyses  of  these 
implants  are  shown  in  Fig.  8. 

These  calibration  lines  are  even  more  valuable  when  a  quantitative  SIMS 
analysis  in  either  a  single  matrix  or  a  layered  multimatrix  sample  is 
desired.  Because  these  types  of  calibration  lines  are  highly  reproducible, 
dopant  distributions  in  a  homogeneous  group  IIIA-VA  compound  matrix  can  be 
quantified  using  the  calibration  lines  and  a  single  standard  prepared  from  any 
of  the  group  IIIA-VA  compound  matrices.  There  is  no  need  to  make  a  separate 
standard  for  each  matrix.  As  shown  previously  for  AlxGaj_xAs  matrices,  these 
calibration  lines  also  make  possible  the  quantitative  SIMS  analysis  of  layered 
nultimatrix  samples.  When  the  criteria  necessary  for  SLIC  are  satisfied  (for 
a  homologous  series  of  compounds  such  as  AlxGa^_xAs,  GaAsxPj_x,  and  InAsx 
Sbi_x),  very  complex  layered  multimatrix  samples  may  be  quantitatively 
analyzed.  When  SLIC  is  not  applicable,  simpler  structures  with  just  a  few 
well  defined  interfaces  may  be  quantitatively  analyzed  by  treating  the  inter¬ 
faces  as  linear  concentration  gradients  from  one  matrix  (of  predetermined 

composition)  to  another  (the  matrix-by-matrix  approach).  Both  approaches  have 

Q  + 

been  applied  to  the  hypothetical  MBE  structure  illustrated  in  Fig.  9a.  A  Be 
concentration  plateau  approximately  0.1  urn  wide  and  1-2  x  10*®  atom/cm^  high 
was  grown  by  MBE  while  the  matrix  was  linearly  changed  from  GaAs  to 

Q 

Al.30Ga.70As.  In  the  uncorrected  profile.  Fig.  9b,  the  Be  distribution 
(dashed  line)  resembles  a  sharp  spike  rather  than  a  plateau.  In  addition,  the 
thickness  and  maximum  concentration  of  the  plateau  cannot  be  determined.  The 
matrix  corrected  depth  profiles  using  the  matrix-by-matrix  approach  and  SLIC 
are  presented  in  Figs.  9c  and  9d,  respectively.  The  shape  of  the  plateau 
differs  slightly  between  the  two  versions  of  the  corrected  profile.  However, 
both  types  of  corrections  permit  a  dramatic  improvement  over  the  uncorrected 
profile  in  the  determination  of  plateau  thickness  and  maximum  concentration. 


L_ 


TIME  (arbitrary  units) 


TIME  (arbitrary  units) 


SIMS  Analysis  of  UB+,  24Mg+,  and  7Li+  Implants  Through  a 

CaAs/AlQ^jGaQ^yAs  Structure,  (a)  Depth  profiles  of  ^B"1”  (...) 
an(I  As  (”■■”)  (interface  at  31  time  units);  (b)  depth  profiles 

Mg*  ( - ) »  +  ( - )  and  7^As+  (...)  (interface  at  20  time 

units).41 
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CONCENTRATION 


Fig.  9.  SIMS  Analysis  of  a  Be  Plateau  at  a  Graded  GaAs/Alg^Cag^ ^As  Inter¬ 
face.  (a)  Hypothetical  structure  Be  (■—”)  (1.2  *  10  atom/cm  peak 
cone.)  and  A1  (— — )  (6.7  x  10Z1  aton/cm3  peak  cone.). 
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Fig.  9.  SIMS  Analysis  of  a  Be  Plateau  at  a  Graded  GaAs/Alg^Gag^  7AS  Inter 
face,  (b)  Uncorrected  profile  of  Be+  ( - )  and  ^As+  (— — ). 
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Fig.  9.  SIMS  Analysis  of  a  Be  Placeau  aC  a  Graded  GsAs/AIq^ j^a0. 7^  *nter“ 

face,  (c)  Corrected  profiles  of  Be  (— )  (2  xlO*®  atom/cm  full 

22  3 

scale)  and  A1  (— )  (1  xlO“  atoo/cm■,  full  scale)  obtained  assuming 
a  linear  concentration  gradient  at  the  interface. 
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CONCENTRATION 


(d) 


Fig.  9. 


SIMS  Analysis  of  a  Be  Plateau  at  a  Graded  GaAs/Alg^ 3GaQ^ yAs  Inter¬ 
face.  (d)  Corrected  profiles  of  Be  (— )  (2  xlO*®  atom/cm^  full 
scale)  and  A1  (— )  (1  xlO^  atom/ca^  full  scale)  obtained  using  the 
program  SLIC. 
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VII.  FUTURE  TRENDS 


These  Initial  investigations  into  the  calibration  of  the  SIMS  matrix 
effects  for  the  quantitative  analysis  of  multilayers  of  IIIA-VA  semiconductors 
have  established  the  groundwork  for  future  investigations.  In  particular, 
five  important  points  have  been  made.  First,  both  ion  yields  and  sputtering 
yields  were  shown  to  vary  linearly  with  the  matrix  composition  of  AlxGa^_xAs 
and  similar  series  of  IIIA-VA  semiconductors.  Second,  matrix  calibration 
lines  were  shown  to  have  superior  precision  and  reproducibility  when  relative 
ion  yields  and  relative  sputtering  yields  were  used  instead  of  useful  ion 
yields,  relative  sensitivity  factors,  and  sputtering  yields.  Third,  a 
procedure  for  the  quantitative  SIMS  analysis  of  multilayer-multimatrix  samples 
has  been  established  and  incorporated  into  the  program  SLIC.  Fourth,  the 
investigation  into  the  factors  influencing  matrix  effects  in  group  IIIA-VA 
semiconductors  revealed  that  relative  ion  yield  variations  with  matrix  are  a 
linear  function  of  the  matrix  composition  and  are  influenced  primarily  by  the 
oxygen  affinity  of  the  matrix  constituents.  Fifth,  for  elements  in  the  same 
column  of  the  periodic  table,  the  influence  of  matrix  effects  on  relative  ion 
yields  was  shown  to  increase  as  the  first  ionization  potential  of  the  analyte 
Increased. 

Although  these  accomplishments  represent  a  significant  advance  in  the 
development  of  quantitative  SIMS  analysis,  matrix  effect  calibration,  depth 
profile  correction,  and  the  mechanisms  of  matrix  effects  should  all  be 
examined  in  much  greater  detail.  Because  of  the  large  inconsistency  in  the 
Ion  and  sputtering  yields  obtained  by  any  given  laboratory,  these  parameters 
cannot  be  evaluated  critically  in  the  open  literature.  As  a  result,  progress 
in  the  understanding  of  the  sputtering  and  ion  emission  processes  has  been 
extremely  slow.  This  difficulty  could  be  minimized  through  the  use  of  the 
more  reproducible  relative  ion  yield  and  sputtering  yield  values.  Thus,  it 
should  be  common  practice  to  report  Rr  and  RS  values  for  all  materials  being 
examined  by  SIMS.  When  so  doing,  a  single  material  (perhaps  silicon)  should 
be  established  internationally  as  the  standard  matrix  for  normalization. 


Thus,  ion  yields  and  sputtering  yields  obtained  by  different  laboratories 
could  be  more  readily  compared  and  correlated.  In  particular,  the  Rrs  and  RSs 
of  homologous  matrices,  such  as  metal  silicides  and  germanides,  should  be 
examined  for  linear  relationships  like  those  established  for  the  IIIA-VA 
semiconductors. 

Any  procedure  for  quantifying  SIMS  depth  profiles  through  multimatrix 
films  must  consist  of  three  steps:  matrix  identification,  sputtering  rate 
determination,  and  ion  yield  calibration.  Although  each  of  these  steps  can  be 
performed  simultaneously  by  SLIC,  the  procedure  used  in  SLIC  to  determine  the 
matrix  composition  requires  the  presence  of  a  single  element  (such  as  As  in 
AlxGa1_xAs)  of  known  concentration  in  each  of  the  layers.  This  requirement  is 
not  generally  satisfied.  Thus,  a  general  method  of  matrix  determination 
through  layers  is  required  to  complement  the  matrix  effect  calibration.  One 
possible  solution  is  the  use  of  a  second  analytical  technique  (such  as  RBS  or 
AES)  to  determine  matrix  composition  and  structure  while  matrix  calibrated 
SIMS  is  used  to  determine  the  dopant  distributions.  As  previously  mentioned, 
small  mismatches  in  the  matrix  and  trace  element  depth  distributions  can 
result  in  large  errors  in  the  matrix  corrected  trace  element  distributions. 
Thus,  the  critical  factor  in  such  a  dual  analysis  approach  is  a  coherent  depth 
distribution  match  between  the  two  data  sets.  Methods  of  ensuring  an  adequate 
registration  between  the  two  data  sets  need  to  be  developed. 

The  reproducible  nature  of  Rrs  should  be  exploited  to  determine 
mechanisms  for  the  influence  of  matrix  effects  on  practical  ion  yields.  The 
data  Indicate  that  t's  are  a  function  of  matrix  composition  and  also  the 
surface  concentration  of  reactive  primary  species  (oxygen  and  cesium).  The 
surface  concentration  of  reactive  species  depends  on  the  availability  of  these 
species  through  surface  adsorption  and  implantation  during  ion  bombardment, 
and  also  on  the  oxygen  affinity  of  the  matrix.  However,  little  work  has  been 
done  to  determine  the  relative  importance  of  these  factors  and  how  experi¬ 
mental  parameters  may  be  altered  to  permit  matrix  calibration  for  a  wider 
variety  of  samples.  The  contribution  of  adsorption,  implantation,  and  matrix 
composition  must  be  deconvoluted  so  that  matrix  effects  observed  for  a  variety 
of  samples  and  analysis  conditions  may  be  understood. 
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LABORATORY  OPE RATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  allltary  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experiaental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  Investigations  is  the  technical  staff's  wide-ranging  expertise  and  Its 
ability  to  stay  current  with  new  developments.  This  expertise  Is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  Individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  including  chemical  kinetics, 
spectroscopy ,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  sppllcatlons ,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solld-stste  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  mlcrowave/allllmeter  wave  measurements, 
diagnostics  and  radlometry,  mlcrowave/mllllmecer  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-parttcle  Interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics.  Infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  Ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
Instrumentation. 
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